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ABSTRACT 

The anisotropic nature of solar wind magnetic turbulence fluctuations is investigated scale-by-scale using high 
cadence in-situ magnetic field measurements from the Cluster and ACE spacecraft missions. The data span five 
decades in scales from the inertial range to the electron Larmor radius. In contrast to the inertial range, there 
is a successive increase towards isotropy between parallel and transverse power at scales below the ion Larmor 
radius, with isotropy being achieved at the electron Larmor radius. In the context of wave-mediated theories 
of turbulence, we show that this enhancement in magnetic fluctuations parallel to the local mean background 
field is qualitatively consistent with the magnetic compressibility signature of kinetic Alfven wave solutions of 
the linearized Vlasov equation. More generally, we discuss how these results may arise naturally due to the 
prominent role of the Hall term at sub-ion Larmor scales. Furthermore, computing higher-order statistics, we 
show that the full statistical signature of the fluctuations at scales below the ion Larmor radius is that of a single 
isotropic globally scale-invariant process distinct from the anisotropic statistics of the inertial range. 
Keywords: plasmas — magnetic fields — solar wind — turbulence — data analysis 



1. INTRODUCTION 

In-situ measurements of fields and particles in the interplan- 
etary solar wind provide unique observations for the study of 
a turbulent plasma (Tu & Marsch 1995; Horbury et al. 2005; 
[Bruno &~Carbone 2005) . As in neutral fluid turbulence, an 
inertial range of magnetohydrodynamic (MHD) turbulence 
is indicated by the observations of a power spectral density 
(PSD) exhibiting a power-law form over a large range of 
scales. This power-law is a manifestation of a scale-invariant 
turbulent cascade of energy from large to small scales. Unlike 
neutral fluids where dissipation is carried out by viscosity, in- 
terplanetary space plasmas are virtually collisionless. On spa- 
tial scales characteristic of the ions there is a transition in the 
PSD from the inertial range at lower frequencies to a steeper 
power-law spanning up to two decades in scale to the elec- 
tron scales (Sahraoui et al. 2009, Kiyani et al. 2009b). This 
second interval of scaling in the PSD was dubbed the dissi- 



pation range in a nalogy with hydrodynamic turbulence ( Lea- 
|mon et al.|1999| . The nature of the fluctuations on these ki- 
netic scales, the mechanisms by which the turbulent energy 
is cascaded and dissipate d, and the p ossible role of dispersive 
linear wave modes ( Leamon et al.|1999; Bale et al. 2005 ; Gary] 
[et al. 2008, 2010| are all hotly debated as vital ingredi ents of 
any future modelof this dissipation range ([Sc hekochihi n et al. 



2009; Ale xandrova et al.|2008[ [Howes et al. 2008; Sa hraoui 
et al.|2010i 



Anisotropy, with respect to the background magnetic field 
is a central feature of plasma turbulence in the solar wind 
IMatthaeu s et al.|1996[JOsman & Horbury|2007[ [Narita et al.| 



dM 
|2D 



10; Sahr aoui et al.|26lO|l , with the fluctuating components 
transverse and parallel to the background magnetic field dis- 
playing manifest differences in bo th their dynamics and statis- 
tics. The seminal study of [Belcher & Davis ( 1971| l used 
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Mariner 5 observations to investigate the field variance ten- 
sor projected onto an orthonormal 'field-velocity' coordinate 
system. They found that the majority of the power is in the 
transverse fluctuations with a ratio between transverse and 
parallel (compressible) components ~ 9 : 1. As first seen in 
high cadence WIND spacecraft data (Leam on et al.|[T998b[ ), 
this ratio decreases to ^ 5 : 1 in the dissipation range, in- 
dicating enhanced magnetic compressible fluctuations in this 
sub-ion Larmor scale range. This result showing an enhanced 
level of magnetic compressibility was repeated with a large 
set of ACE spacecraft data intervals (Hamilton et al. 2008 1 
and also in much higher cadence data from the Cluster space- 
craft (Alexandrova et al. 2008). 

In this article we conduct a scale-by-scale study of the 
anisotropy in transverse and parallel fluctuations over five 
decades in temporal scales from a few hours to tens of mil- 
liseconds. Our background magnetic field is scale-dependent, 
computed self-consistently with the scale-dependent fluctua- 
tions as in ( Horbury et al.^2008, Podesta. 2009) i.e. a local - 
as opposed to globa l - background magnetic field. We show 
that th e findings of ( Leamon et al.| 1998b Alexandrova et al. 
|2008| l are in fact the result of a successive scale-invariant re- 
duction in the power ratio between the two components as we 
move to smaller scales in the dissipation range. This reduction 
results in isotropy between all three components of magnetic 
field fluctuations (two transverse and one parallel) at the elec- 
tron Larmor radius p,, i.e. equipartition of magnetic energy 
between all three magnetic field components. Using a Hall 
MHD model we provide a simple description of how the Hall 
term in the generalised Ohm's law is responsible for the rise 
of parallel (compressible) magnetic fluctuations and the resul- 
tant isotropy at scales around p,,. In addition, by comparing 
with linear kinetic solutions of the V lasov equ ation and com- 
puting the magnetic compressibility ( [Gary & Smith 2009 ) we 
show that this scale-invariant reduction in the power ratio is 



2 



also qualitatively consistent with the transition from shear or 
highly oblique Alfvenic fluctuations in the inertial range pos- 
sessing no magnetic compressibility, to kinetic Alfven wave- 
like fluctuations which rapidly develop magnetic compress- 
ible co mponents at scales around the ion Larmor radius ( jGary| 



& Smith 2009 ; Sahraoui et al. 2012 ). This nature of the mag- 
netic compressibility is also suggested in the recent work by 



Salem et al. ( 2012[ l who, using a global background field, 
also computed the magnetic compressibility and compared it 
with wave solutions of the linearised warm two-fluid equa- 
tions. Importantly, by calculating higher order statistics and 
the probability density functions (PDF) we show for the first 
time that the dissipation range magnetic fluctuations are char- 
acterised statistically by a single isotropic signature; in stark 
contrast to the anisotropic inertial range. 

We stress that the anisotropy discussed in this article is in 
the full vector direction of field fluctuations, (5B, projected 
parallel and transverse to the local magnetic field directio n 
(also known as variance anisotropy ( [Belcher & Davis|1971| l). 
This is in contrast to the anisotropy measured in the vector 
r which characterises the spatial scale of the fluctuations (or 
wavevector k in Fourier space). In the context of the latter 
anisotropy, the spatial variations in the magnetic field that we 
measure are a sum of all the components of k which are pro- 
jected along the spacecraft trajectory through the solar wind 
i.e. along the solar wind velocity unit vector -Vsw; these are 
then Doppler shifted (via the Taylor frozen-in-flow hypothe- 
sis) to spacecraft frame frequencies, /, resulting in a reduced 
spectrum P(f, Ogv) = J c/-'kP(k)(5(27r/-k • Vsw), where the S- 
function is doing the 'reducing', and Obv shows the explicit 
dependence of the reduced spectra on the angle between the 
background magnetic field and th e bu lk solar wind velocity 
vectors (see |Forman et al.] ( |2011[ ) and [Fredricks & Coroniti| 
([T976j) for further details). The full k spectrum, /"(k), can only 
be measured by multi-spacecraft techniques which take into 
account th e full 3D spatial variation of the magnetic field e.g. 
^-filtering ( jPingon & Lefeuvre|199r) . In the context of the re- 
duced spectrum, the magnetic field observations investigated 
in this article are dominated by wavevector components which 
are strongly oblique (70° - 90°) to the background magnetic 
field; as indicated by the angle between the background mag- 
netic field and solar wind velocity vectors. In such intervals, 
the inertial range is ubiquitously seen to have a spectral in- 
dex of ~ -5/3, whilst in the dissipation range this steepens to 
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-2 .8 (|Kiyani et al |[2009 b; Alexandrova et al. 2009; Chen] 
"aL]|20I0[ [Sahraoui et aLpOlO) ). Also, for such int ervals; 



multi s pacecraft observation s of the full f (k) (N arita et al.[ 
[2010[ Sahraoui et al.[[20T0| ) suggest that most of the power 
IS m the k± rather than the k^<^ components, suggesting that 
highly oblique angled components of k are very much repre- 
sentative of most of the magnetic field fluctuation power. 

A brief synopsis of the paper is as follows: Section |2] de- 
scribes the spacecraft data and summarises the plasma param- 
eters involved, as well as the analysis methods used; Section[3] 
forms the bulk of the paper and is split between a.) a study and 
discussion of the rising magnetic compressibility and compo- 
nent anisotropy in the dissipation range and b.) an anisotropic 
study of the higher-order statistics. Finally Section |4] sum- 
marises our findings and brings together the insights obtained 
from these results, concluding with an outlook for future in- 
vestigations. Details of the Undecimated Discrete Wavelet 
Transform that we use to decompose the magnetic field into 
scale-dependent background field and fluctuations, for use in 
this study, are described in Appendix [A| 



2. OBSERVATIONS AND METHODS 

We discuss an interval of quiet ambient solar wind for 
which there are observations from both the Cluster and ACE 



spacecraft missions (Escoubet et al. 1997 ' ^Garrard et al. 1998| l 
(all in GSE coordinates). ITie Cluster ( spacecraft 4) interva l 
at 450 Hz cadence (same interval as in (Kiya ni et al.|2009"b i) 
is of an hour duration 2007/01/30 00:10-01:10 UT when the 
instruments were operating at burst mode, and will primarily 
be used to study the dissipation range at spacecraft frequen- 
cies above 1 Hz. We construct a combined data set from the 
DC magnetic field (sampled at 67 Hz) of the Flux Gate Mag- 
netometer (FGM) for frequencies below 1 Hz, with the high 
frequency (oversampled at 450 Hz) search-coil magnetome- 
ter data from the STAFF-SCM experiment for frequencies 
above 1 Hz, using the wavelet reconstruction method ( Alexai>] 
drova et al.[2004, ,Chen et al.,2010) . The ACE magnetometer 



(MAG) interval at 1 Hz cadence is over two days 2007/01/30 
00:00 UT to 2007/01/31 23:59 UT and samples the inertial 
range at spacecraft frequencies below 1 Hz. The ACE in- 
terval was only needed in order to obtain better estimates of 
the higher-order statistics in the inertial range (as the Cluster 
FGM interval is relatively short for this), and thus will only 
be restricted to this study in the penultimate section of the pa- 
per; for the other studies of the inertial range, the FGM data 
interval was used. 

Both these intervals are in stationary fast wind (~ 
667 km i"') with similar plasma parameters and are suffi- 
ciently large data sets so that the sample size variance er- 



rors in the computed statistics are negligible (Kiyani et al. 
'2009a). The following plasma quantities are from the Clus- 
ter FGM, CIS, PEACE and WHISPER instruments: average 
magnetic field B ~ 4.5 nT, electron plasma /3g ~ 1.2, plasma 
density n^ ~ 4 cm~^, Alfven speed Va — 50 km s~\ perpen- 
dicular ion temperature Ti± ~ 24 eV, electron temperature 
~ 22 eV, and ion and electron Larmor radii p, ~ 110 km 
and ~ 1 .7 km respectively. The interval is free from any ion 
or electron foreshock effects at the Earth's bowshock. More- 
over, using high-resolution proton moments computed from 
the 3DP instrument ( [Lin et al. 1995) onboard the Wind space- 
craft at lAU and at the same interval as Cluster, the proton 
plasma parameters (parallel and perpendicular proton temper- 
atures Tp\\ ~ 26 eV and Tp± ~ 34 eV, parallel proton plasma 
~ 1 .2) indicate that the intervals are stable to proto n pres- 
sure anisotropy-driven instabilities ( Bale et aL][2009| l. From 
flie ACE S WEPAM instmment the ion plasma beta /3, ~ 1 .5 

showing that the ion-inertial le ngth A,- c± p,-. 

There is a growing opinion (Chapman & Hnatl2007 Hor- 
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a local scale-dependent mean field consistent witH the scale 
dependent fluctuations, rather than a large scale global field, 
which should be used in studies of anisotropic plasma tur- 
bulence. This is not only self-consistent but also assures 
that there are no significantly large spectral gaps between 
the frequencies of the fluctuations being studied and the 
mean-fields being projected upon. In accordance with this 
approach we use the Undecimated Discrete Wavelet Trans- 
form (UDWT) method to decompose the fields into local 
scale-dependent background magnetic fields and fluctuations, 
B(f ,/) and (5B(f ,/) respectively; where / explicitly shows the 
frequency or scale dependence - the parallel and transverse 
fluctuations, ^B(||/j^)(f,/), are then obtained from these. The 
background to this particular wavelet method and brief details 
of the algorithms are described in Appendix [A| 
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Figure 1. (Upper panel) PSD (from Cluster) of the transverse and parallel 
components spanning the inertial and dissipation ranges. Standardised sam- 
ple size variance errors are smaller than the markers. The search-coil sensitiv- 
ity floor PSD is obtained from the z-component (spacecraft SR2 coordinates) 
of a very quiet period in the magnetotail lobes (2007/06/30 15:00-15:05 UT) 
as a proxy for the instrumentation noise. (Lower panel) Ratio of parallel over 
transverse PSD. Horizontal dot-dashed line indicates a ratio of 1 /3 where 
isotropy in power occurs. Vertical dashed and dashed-dotted lines indicate p, 
and pe respectively, Doppler-shifted to spacecraft frequency using the Taylor 
hypothesis. The litde blip in PSD\~^ at ~ 0.25 Hz (and in PSDJPSD±_) is 
due to the residual spacecraft (~ 4 seconds) spin-tone in the FGM signal; it 
is more noticeable in PSZ)|| due to the lower power in pai'allel fluctuations in 
the inertial range. 

3. RESULTS AND DISCUSSIONS 

3.1. Power isotropy and enhanced magnetic compressibility 

In keeping with Parseval's theorem for the conservation of 
the L^-norm (energy conservation) the wavelet PSD for the 
parallel and transverse magnetic field components is given by 
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where SBj_(tjJ) = J 5B\^{tj,f) + 5B\^{tjJ) is the magni- 
tude of the total transverse fluctuations at time tj and fre- 
quency /, A is the sampling period between each mea- 
surement, and is the sample size at each frequency /. 
For the Cluster interval the PSD\\ and PSD± are shown in 
fig[T] The spectral indices obtained are ~ -1.62 ±0.01 and 
~ -1.59 ±0.01 for parallel and transverse components re- 
spectively, in the inertial range; and ~ -2.67 ± 0.01 and 
~ -2.94±0.01 for parallel and transverse components respec- 
tively, in the dissipation range. The lower panel of fig|T|shows 
that not only do these results r ecover the ^9:1 amsotropy 
ratio of [Belcher & Davis] ( | 1971, ), they also show that the de- 



crease in the anisotropy observed by (Leamon et al.|| 1998b 
[Alexandrova et al. [[20081 [Hamilton et al.||2008p in the dissi- 
pation range is actually a scale-free progression to isotropy. 
This progression of the anisotropy in the power ratio PSDy± 
(a measure of magnetic compressible fluctuations) begins at 
the spectral break (spacecraft frequency ^ 0.25 Hz), just be- 
fore the calculated p,, and follows the power-law relationship 
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to Pe, where isotropy in power between 



the three components (parallel and two transverse compo- 
nents) is achieved. This isotropy corresponds to a value of 
PSD^^ /PSD± = 1/2 and is indicated in the lower panel of fig. 
[T| Although this enhancement of parallel, or compressible, 
mictuations in the dissipation range ha s already been com 
mented upon by various authors (Leamon et al. 1998b Hamil 



ton et al. 2008 Alexandrova et al.|2008jl, it is n ormaffy shown 
to be nearly constant (apart from in ( [Salem e t al. 2012)). To 
our knowledge this is the first time that an observation of 
isotropy has been noted to occur at kp^ ~ 1; although it is 



also strongly suggested in the PSDs in Sahraoui et al. (20101. 



The lower panel of fig. [T] is calculated from the ratio of the 
averages of the parallel and transverse fluctuations to show a 
measure of the anisotropy. As such it does not constitute a 
proper estimate of the ensemble average of the anisotropy. A 
proper measure would be to take the ratio of 5Bl(tj,f) and 
5B\(tj,f) at each time fy, and then average over this ensem- 
ble of realisations of the anisotropy. However, this is prone to 
large errors induced by very large spikes caused by purely, or 
nearly pure, parallel fluctuations which result in divisions by 
a very small number. To overcome this problem and to obtain 
a proper ensemble averaged anisotropy measure, we compute 
an alternative and robust metric of the magnetic compressibil- 



ity (similar in form to the expressions in (Gary & Smith 2009 
[Alexandrova et al.|2008| l) defined as 
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i.e. the compressibility is calculated locally from the time- 
dependent fluctuations and then averaged. Converting the 
spacecraft frequency into wavenumber using the Taylor hy- 
pothesis and normalising by the averaged ion gyro radius for 
the interval, figure l2| shows C\\{kpi) computed using the lo- 
cal scale-dependentmean field. We have binned the latter 
into 10° angle bins (angle between y sk and e||(fy,/)) to show 
which components of the k variation we are measuring with 
respect to the local background (scale-dependent) magnetic 
field. On the same plot, and for comparison, we have also 
included the calculation of the magnetic compressibility us- 
ing a global background magnetic field - the latter consists of 
the mean average of the magnetic field vector over the whole 
interval being studied. The angle of the solar wind veloc- 
ity to the background magnetic field using this global mean 
field would correspond to wavevector components at an an- 
gle of ~ 75° to the background field, both within the inertial 
and dissipation ranges. From these plots of C\\ {kpi) we can 
see that there is not only a large difference between C\\ calcu- 
lated using a local as opposed to a global field, there is also 
no significant difference between C\\ calculated using a local 
field for the separate k-component angle bins. Also, simi- 
lar to the lower panel of fig. [T] in all the curves we again 
see the enhancement of the magnetic compressibility - with 
isotropy between all three components of the vector magnetic 
field fluctuations being reached at kpe ~ 1 . 
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Figure 2. Magnetic Compressibility C| | (from Cluster), computed as in eqn. 
(21, against the normalised wavenumber for both local (scale-dependent) and 
global (computed over the whole interval) background magnetic fields. The 
normalised wavenumber was computed assuming the Taylor Hypothesis. Nu- 
merical solutions of the linearised Vlasov equation, corresponding to the ki- 
netic Alfven wave mode, (computed using the solar wind parameters for this 
interval) are also shown in comparison to illustrate the rising magnetic com- 
pressibility possessed by these modes at the ion-gyroscale. The horizontal 
dashed line at C| | = 1/3 indicates the level at which full isotropy of the mag- 
netic field fluctuations occur. As in fig. ifTJ the statistical sample size variance 
errors are smaller than the markers, anothe effect of the spacecraft spin is 
clearly seen at kpt ~ 0.25. 

Lastly, one should note from fig. |2]that the difference be- 
tween the magnetic compressibility computed with global and 
local scale-dependent background magnetic fields decreases 
as we move to smaller scales. At a first glance this is a non- 
intuitive result, as we would expect that the local field increas- 
ingly mimics the global field at larger scales. However, this 
can be easily resolved by the fact that if the fluctuations are 
isotropically distributed, it does not matter which basis one 
is projecting the fluctuations in - they will look the same in 
all bases. So if the fluctuations are becoming more isotropic 
towards the smaller scales, the less will be the difference be- 
tween the global and local background magnetic field bases. 

3.1.1. Rising magnetic compressibility in kinetic Alfven wave 
solutions of the linearised Vlasov equation 

Theories of plasma turbulence which advocate that a turbu- 
lent energy cascade within the dissipation ran ge is mediated 
by the various linear wave modes of a plasma ("Leamon et al. 
T999||Schekochihin e t al. 2009 ; Gary et al. 2008 ), suggest pre- 
dictions for C|| (Gary & Smith|2009||Sahraoui et al.|2012| . As 
the wave modes advocated to cascade energy at tnese scales 
are dispersive in nature (kinetic Alfven and Whistler waves), 
many of the proponents of such theories suggest that what has 
been dubbed the dissipation range for purely histor ical rea- 
sons, should actually be called th e dispersive range ( jStawicki 
let al.|200H|Sahraoui et al.|2012| l. Within the context of sub- 
ion Larmor scales, there is a growing body of work showing 
that the fluctuations at these sc ales share the characteristic of 
kinetic Alfven wave s (Leamo n et al.|199 8b, 1999; Bale et aL 
^005] [Howes et al.||2008, .Sahraoui et al..,2009, .2010, ,Salem 



waves into kinet ic Alfven waves, with a strong compressive 
component (HoUweg 1999 1 and propagating at highly oblique 
(near-perpendicular) angles to the background magnetic field. 

In figurel2|we overlay the observational results of C\\ by pre- 
dictions orU||(A:p,) from numerical solutions of the linearized 
Vlasov-equation for the kinetic Alfven wave mode at a highly 
oblique (85°) angle and at a virtually perpendicular (89.99°) 
angle, using the same bulk plasma field and particle parame- 
ters of our data. These numerical solu tions were plotted using 
the WHAMP code ( Ronnmark l 1 982) l which assumes a strong 
guide field i.e. a global background magnetic field. The plot 
shows that our results are qualitatively consistent with C\\ pre- 
dicted for kinetic Alfven waves. This enhancement of the 
compressive component in the kinetic Alfven wave is due 
to the coupling of the A lfven mode to the pu rely compres- 
sive ion acoustic mode ( [Leamon et al.|[T999] l. In the iner- 
tial ran ge below kpi ~ 1 both of the se modes are decoupled 
( [Howes et al.„2()0^ |Schekochihin et al. 200 9). In the con- 
text of such theories, the 10% compressible component that 
we see in our results in the inertial range is most likely due 
to the presence of magnetosonic compressible modes in the 
data, all at highly oblique (nearly perpendicular) angles to the 
background magnetic field. The modes can be in the form 
of slow modes, fast modes or non-propagating (a; = 0) pres- 
sure balanced structures advected by the flow (mirror modes), 
all of which have purely parallel (compressible) fluctuations 



et al.|2012[ i. Indeed, the enhancement of the magnetic com 
pressibility seen above is strongly consistent with the transi 
tion, around kpi ~ 1, of purely incompressible shear Alfven 



(ISahraoui e t al.[ 2012). Although a recent paper by Howes 
et al. (201Ia| shows that slow modes or pressure balanced 
structures constitute the dominant part of the compressible 
fluctuations in the inertial range, fast magnetosonic modes 
may also exist. This can be supported by the fact that these 
fast modes are shown to be considerably undamped (com- 
pared to slow modes) in the linear theory at highly oblique 
angles of propagation and high plasma beta (/?,■ > 1) (Sahraoui] 
et al. 2012), conditions very relevant to the Solar Wind. In fig- 
ure[2] the departure of our results from the theoretical solution 
at kpi < 1 can thus possibly be explained by the fact that the 
theoretical curve is only for kinetic Alfven waves - an addi- 
tion of the above purely compressible magnetosonic modes 
might be able to reduce the difference significantly. In addi- 
tion the theoretical plot is calculated assuming a strong global 
background magnetic field. Importantly, the theoretical plot 
for kinetic Alfven waves shows that the magnetic compress- 
ibility should start to ris e before kpi ^ 1 whi ch is also what our 
results show. Recently, Salem et al. p012| l showed a similar 
result for magnetic compressibility (their definition of C|| not 
being a squared quantity as in this article) with better agree- 
ment with the theoretical C|| curve for kinetic Alfven waves 
than the res ults presented here. The differences in our results 
and those of [Salem et al. ( [2012| l could be due to the intervals 
analysed: our interval is fast wind with /3, ^ 1 .5, whereas their 
interval is slow wind with /3, ^ 0.5; that they used two-fluid 
linear warm plasma relations, instead of the linear Vlasov so- 
lutions used here; or that their component fluctuations were 
projected onto a global background magnetic field instead of 
a local scale-dependent one. To put both these results in their 
proper place, requires a large data survey with intervals that 
sample more broad and comprehensive plasma conditions. 

3. 1.2. The role of the Hall term in the enhancement of magnetic 
compressibility and component isotropy 

In the more general non-linear setting, this enhancement of 
magnetic compressible fluctuations can be explained by the 
increased prominence of the Hall term, j x B, in the gener- 
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alised Ohm's Law, where j is the current density arising from 
unequal ion and electron fluid velocities. To show this in the 
dynamical evolution of the magnetic field, we will turn to the 
Hall MHD model, as it is the si mplest physical m odel of a 
plasma with unmagnetised ions ( [Shaitdi & Sh ukla 2009). Al- 
though Hall MHD might not be an accurate representation of 
the full kinetic physics, especially the linear kinetic physics 
in hot plasmas ( Ito et al.]|2004 Howes|[2009) , it serves as a 
simpler nonlinear model to explain the nsmg magnetic com- 
pressibility and the isotropy that we observe here. Primarily, 
we will look to the Hall-MHD induction equation (in SI units) 

^ = Vx(vxB)+— Vx(B-VB) + ?yV-B, (3) 
at iierig 

where v is the bulk (ion) velocity, rj is the magnetic diffusiv- 
ity, fi is the permeability, e is the electron charge and is the 
electron number density. The first term on the RHS of equa- 
tion Q is the convective/dynamo term, the second term is the 
Hall term and the last term is the diffusive term. We have ne- 
glected the effects of any electron pressure gradients here to 
simplify the discussion, and as these are OiSn^/n^) their ef- 
fects are considered small. We have kept the diffusive term 
- which is normally negligible for these solar wind parame- 
ters - in order to retain an energy sink in the equations, but 
also for the possibility of a renormalized, or turbulent, mag- 
netic diffusivity; however, we will not be making use of this 
term in the following arguments. Although other terms, such 
as the electron inertia that are very relevant to high frequency 
modes e.g. Whistler modes, become relevant at scales kX^ ^ 1 
(where A^. is the electron inertial length), we stick to the Hall 
MHD model above and focus on the effects of the Hall term 
with respect to the convective term. Importantly, in (3 ^ C(l) 
plasmas (A/ ^ p, as A,- = pij \fP), such as the solar wind at 
1 AU, as one approaches the ion-Larmor scale p,, the Hall term 
becomes of the order of the convective term (Goossens 2003jjl. 
In the inertial range, at scales above the Larmor radius, it is 
the convective term which dominates and the magnetic field 
is frozen to the ion flow, such that any flow of the plasma per- 
pendicular to the magnetic field affects the evolution of the 
magnetic field. Thus, if we assume that much of the power in 
inertial range magnetic field fluctuations is locally generated 
from the velocity fluctuations (via the non-linear evolution of 
the ion momentum equation) in the form of a dynamo effect 
(V X (v X B) term in ([3])), and is not being passively advected 
in the solar wind from the Sun, then it is natural that most 
of the power in the inertial range will be in transverse fluctu- 
ations - also considered a signature of Alfven ic fluctuations 
( [Belcher & Davis|197I| JHorbury et"^[2005] l. Any residual 
parallel (compressible) fluctuations arise from fluctuations in 
the plasma density (which also contribute to the transverse 
fluctuations). 

If we now approach scales close to the Hall term starts 
to show its effects. This manifests itself in changing the di- 
rection of the fluctuations such that if transverse fluctuations 
dominate, the Hall term will generate parallel fluctuations and 
thus result in an increase in the magnetic compressibility C|| . 
To see this we will focus on the Hall term and assume an over- 
simplification that a.) it is entirely transverse fluctuations 8B\^ 
which dominate the turbulent fluctuations coming from the in- 
ertial range and b.) that fluctuations dominate (i.e. perpen- 
dicular gradients Vj^) - both of which are supported by our 
observations, and studies of the full 3D k-spectrum (Sahraouij 
[et al..2010) . This results in the Hall term contribution to the 



induction equation becoming 
9B|| 



dt 



V^i x(B-V,5B^2) + Vi2x(B-V(5B^i), (4) 



Hall 



where we have explicitly separated the two transverse compo- 
nents to make clear the that the two directions are not parallel 
to each other in the vector cross product. Equation (|4]i clearly 
shows the origin of the enhancement in the parallel magnetic 
compressible fluctuations. In A;-space the RHS of equation ([4| 
would be in the form of the following convolution 

y"fl'3jk-B(k-j,f)(kii x,5Bi2a,f) + k±2X^Biia,f)), (5) 

where the wavevector arguments are shown explicitly. More 
importantly, equations ([4|i and ([5]l also suggest the origin of 
the isotropy that we observe amongst the different fluctuation 
components at scales close to pg, where terms similar to ([4]l 
and ([si are dominant and the convective term becomes neg- 
ligible! At such scales, the Hall term will dominate and will 
also convert any parallel fluctuations which are generated (or 
already present) into transverse ones. This will evolve un- 
til a steady state is reached between the various terms cor- 
responding to the vector components of the magnetic field 
fluctuations i.e. isotropy between all the vector components. 
Interestingly, this simple observation of the role of the Hall 
term in the enhancement of magnetic compressibility and the 
resultant isotropy, also suggests some information on the k- 
anisotropy i.e. k\\lk^. If, say, equation (jsjl was dominated by 
k\\ instead of k^, then it is clear that an inertial range domi- 
nated by transverse fluctuations will transition into something 
which is also dominated by transverse fluctuations, with little 
or no parallel fluctuations - this is not what we observe here. 
This indicates that, in our Hall-MHD formalism, power in k^ 
fluctuations is not dominant. Intriguingly, however, this type 
of argument also does not rule out the possibility that both 
k\\ and k^ fl uctuations are significant at sub-Larmor scales 
( Gary et ari2010 ); the only observational study of the full 3D 
J ^(k) spectra at such scales by Sahraoui et al. (2010), seems 
to indicate that this is not the case. This and more detailed 
calculations of the energy transfer between wavevector triads 
of fluctuations will be the subject of a future investigation, so 
we can determine a more precise nature of such steady states. 
Crucially, measurements of the magnetic compressibility, in 
the context of a Hall-MHD model, could possibly constrain 
the measurements of the k vector anisotropy. 

3.1.3. The transition range 

A transition range be tween the inertial an d dissipation 
ranges was suggested by Sahraoui et al. (20101), which was, 
in the context of dispersive waves, cited as a possible sign of 
where magnetic field energy is Landau damped into ion heat- 
ing. This transition range, comprising just under a decade 
of scales around kpi ^ I is also seen in th e reduced spectra 
in ([Kiyani et aL][2009b| [Chen et al.l[20T0l l and is distinctly 
different from the power-law which follows at sm aller scales 
down to kpe ~ 1 (see fig. 1 in ( Kiyani et al. 2009b I and fig. 2 in 
(Sahraoui et al. 2010)). Our arguments above suggest that the 
competition between the convective and Hall terms could also 
explain such a transition range. In this case, and from looking 
at where the compressible fluctuations begin to rise in figfl] 
and fig|2](A:p, ^ 0.1), it seems that the transition range could 
actually span a much larger range of scales. It is important to 
note here that we do not know what the functional form of the 
transition from the convective dominated regime to the Hall 



6 



dominated regime is. All we know is that in the inertial range, 
far from the spectral break, one can neglect the effects of the 
Hall term; well below the ion-Larmor radius and well past 
the spectral break we can neglect the convective term; and at 
kpi ~ 1 (for^j~ 1) these terms are of the order of each other 
( [GoossenSj2003| l. Using the above arguments, and the latter 
observation about the similar strengths of the convective and 
Hall terms at kp, ~ 1 we can make the following observation: 
at kpi ~ 1 the convective term will provide half of the fluctua- 
tion power (50%) which we assume, from the arguments and 
observations above, will be comprised of 10% fluctuations in 
the parallel direction; the Hall term will contribute the other 
half of the fluctuation power (50%), and 33% of these fluc- 
tuations will be in the parallel direction. This means that the 
magnetic compressibility at kpt ~ 1 will be 



Cii ~ 0.5 X 0.33 + 0.5 X 0.1 ~ 0.22 



(6) 



This value is in excellent agreement with the value of C|| at 
kpi ~ 1 extrapolated from the local field curves in fig|2] 

3.2. Higher-order statistics and intermittency 

We next calculate higher-order statistics given by th e struc- 
ture fu nctions (absolute moments of the fluctuations) (Kiyani' 
|et al.| [2009b) for the different components of the magnetic 
field fluctuations with respect to the local scale-dependent 
backgrou nd magnetic field. The m ''' order wavelet structure 
function ( Farge & Schneider|2006| l is given by 



1 

-T 



5B\\(^){tj,T) 



(7) 



where, as detailed in Appendix|Al r = 2' A : / = {0, 1 , 2, 3, . . .} 
is the dyadic time scale parameter related to the central fre- 
quency / used earlier, and A is the sampling period. Scale- 
invariance is indicated by S'"^j_^{t) cx r''*'"^; where C(m) are the 
scaling exponents. The structure functions and corresponding 
scaling exponents C,{m) are shown in fig. [3] for both the iner- 
tial and dissipation ranges using the ACE and Cluster intervals 
respectively. 

Similar to the results by Ki yani et al.| ( |2009b| l, the higher- 
order scaling in the inertial and dissipation ranges are distinct. 
The inertial range shows multi-exponent scaling as evidenced 
by a nonlinea r C,{m) c haracteristic of solar wind turbulence at 
MHD scales ( |Tu & M arsch 1995 ). In contrast, the dissipation 
range is monoscaling i.e. characterised by a linear (^{m) = Hm 
and a single exponent H. Notably both parallel and transverse 
fluctuations in the dissipation range show this different scaling 
behaviour in the inertial and dissipation ranges. However, the 
scaling behaviour between the different components is also 
distinct; with a more pronounced difference shown in the in- 
ertial range. The difference in the dissipation range exponents 
for parallel and transverse fluctuations is simply reflecting the 
different spectral exponents seen earlier in fig[T] 

Before we discuss these higher-order scalmg results, we 
complete the statistical results by finally looking scale-by- 
scale at the individual probability density functions (PDFs) 
for the transverse {e^\ direction) and parallel fluctuations. It 
is necessary to pick one of the transverse directions as the 
combined (magnitude of) transverse fluctuations are positive- 
definite quantities, and thus do not illustrate any symmet- 
ric/asymmetric character of the fluctuations. There is no a 
priori reason, from the symmetry of the physics, for one 
not to expect the fluctuations to be axisymmetrically dis- 
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Figure 3. (a i.) Transverse and parallel wavelet structure functions of order 
1—5 (from the bottom) and (a ii.) resultant scaling exponents for the inertial 
range using the ACE data interval. Structure functions have been vertically 
shifted for clarity, (b i.) and (b ii.) have similar descriptions to (a i.) and 
(a ii.) but are from data in the dissipation range using the Cluster data inter- 
val. These are the anisotropic generlisations of the differences in the scaling 
behaviour be tween the inertial and dissipation ranges first shown by [Kiyani] 
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Figure 4. Rescaled (standardised) transverse and parallel PDFs of fluctua- 
tions in the inertial and dissipation ranges (with Poisonnian error bars). Nor- 
malised histograms with 300 bins each were used to compute the PDFs. Four 
values of t were used in both ranges: r = { 16, 32, 64, 128} seconds for the 
inertial range from ACE, and r = {0.036, 0.071, 0.142, 0.284} seconds for 
the dissipation range from Cluster. 



tributed in the plane perpendicular to the background mag- 
netic field. However, if we look at the separate PSDs for the 
two transverse components of the magnetic field fluctuations 
(not shown here) we will notice that this symmetry is bro- 
ken and the power in the two transverse components is dis- 
tinctly different. When we constructed our scale-dependent 
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orthonormal bases (AppendixjA) it was natural to involve the 
background guide field as it is ubiquitously known to order 
the physics in magnetised plasmas. The other two directions 
perpendicular to this are relatively arbitrary, and we naturally 
chose the stable mean bulk veloci ty field direction, Yg -w, to 
form these in the manner of Belcher & Davis| ( |1971[ l: an 
orthonormal scale-dependent 'field-velocity' coordinate sys- 
tem. However, the solar wind bulk velocity field picks a pre- 
ferred sampling direction (in A:-space) resulting in the reduced 
spectrum mentioned earlier in the Introduction. In introduc- 
ing the velocity field in such a way, the reduced spectrum 
breaks the transverse axisymmetry of the magnetic fie ld fluc- 
tuations and introduces a measurement bias (see also Turnerl 



et al. ( 2011| who use comparisons with MHD turbulence sim- 
ulations, and a model of transverse waves to show the im- 
portance of the spectral slope in this broken symmetry). In 
our scenario where we have taken a constant velocity field 
direction, Vsw, throughout our entire interval (due to the dom- 
inance of the GSE x-component), and the background mag- 
netic field remains mostly highly oblique to Vsw, this bias in 
the reduced spectra remains nearly the same at all times, and 
effectively results in a simple (constant) shift in the (log) fre- 
quency which makes the two transverse PSDs misalign with 
each other - it is essentially a Doppler shift-like sampling ef- 
fect. The magnitude of the transverse vector is not affected 
by such a bias, so all our results above are ambivalent to this 
symmetry breaking. Although the PSD of the two transverse 
components differ in this respect, we can confirm that the 
standardised (rescaled) PDFs for both transverse components 
are near-identical. Thus, the breaking of the transverse ax- 
isymmetry does not affect the results presented in this paper. 
The PDFs for 6B\\ and SB±i are shown in fig. |4l where 



we have used the self-affine scaling operation VsioBia' ) = 
aP{6Bi,T) to rescale (standardise) the fluctuations by their 
standard deviation so to offer a comparison of the functional 
form of the PDFs for different r. We show overlaid rescaled 
PDFs corresponding to four values of t in the dissipation 
range, r= {0.036, 0.071, 0.142, 0.284} seconds (from Clus- 
ter), and in the inertial range, r = {16, 32, 64, 128} seconds 
(from ACE). Figure |4] shows that in the dissipation range 
the same PDF functional form is obtained over the range of 
T for both parallel and transverse fluctuations - suggesting 
that the dissipation range is in this sense 'process isotropic' 
with a (speculatively) single physical process determining the 
dynamics. This supports our arguments and results above 
which show the growing isotropy in the components of the 
magnetic field fluctuations. This is in contrast to the iner- 
tial range where we see that the transverse and parallel PDFs 
are clearly different as we would anticipate from their distinct 
C,{m) shown in fig. [3](a ii). 

In addition to the differences in the scaling behaviour, 
both inertial and dissipation range PDFs are highly non- 
Gaussian indicating high levels of intermittency. In con- 
trast, the rescaled PDFs of the STAFF-SCM instrument noise 
shown in fig. |5]is clearly Gaussian; confirming that our re- 
sults are robust to noise contamination. Furthermore, we 
have analysed other intervals, 2007/20/01 12:00-13: 15UT and 
2007/20/01 13:30-14:10UT, and find broadly consistent re- 
sults with those presented above {H\\ = 0.78, H±^ = 0.95 and 
//|| = 0.8, i/j^ = 0.94 respectively), but with lower signal-to- 
noise ratios. Please note, we consider 'intermittency' in the 
most general sense, to simply reflect the presence of rare, or 
'bursty', large amplitude fluctuations in our signals - some- 
thing which is manifest in our non-Gaussian PDFs in fig. 




[nl/a] 



Figure 5. Rescaled PDF of B- magnetic fluctuations from tiie instrument 
noise proxy of fig. [T] at the same scales r as in fig. l4]for the dissipation 
range. A Maximum Likelihood Estimate clearly shows the Gaussian nature 
of this noise proxy - very different from the highly non-Gaussian statistics 
shown in figR] We can also confirm that the smallest scales in our data which 
coincide witn the noise floor, in terms of power and Gaussian statistics, are 
excluded from the study presented in this ailicle. 



El This is in co ntrast to the popular definition, advocated 
by |Frisch] ( | 1 995] ), which links intermittency to a dependency 
of the kurtosis k{t) = S*(t)/{S^{t))^ (or another higher or- 
der moment ratio) with scale r - a property of multifractal 
(multi-exponent) scaling i.e as seen for the inertial range in 
our results. In this latter definition, the constancy of the kur- 
tosis with scale indicates 'self-similar' signals such as our 
monoscaling signature seen in the dissipation range. There 
is no need to plot the kurtosis for our data as this can be very 
easily seen if we use the scaling relations hip of the m^^-order 
structure function S"'(t) = t<^'"^S"'{1) (see (iKiyani et al.|2006 1 
for notation) in the definition for the kurtosis. This will re- 
sult in k{t) = t'^'*^'^'^^^^ k{1) which, if C(m) is a linear function 
of m (C(«^) = Hm), is independent of r; as is the case for our 
results in the dissipation range. However, even though our re- 
sults show a monoscaling signature in the dissipation range, 
there is still a small nominal dependence on the scale r for 
k(t). This could be simply due to a statistical artifact asso- 
ciated with finite sample size and the fact that the kurtosis is 
sensitive to very large fluctuations. Our very largest events at 
the tails of our distribution are not statistically well-sampled - 
an un-avoidable pit-fall of heavy-tailed distributions, as seen 
by the large errors at the tails of our PDFs in fig|4] 



3.2. 1 . Comparison with other works on scaling at kpi > 1 

We now compare our results for higher-order statistics in 
the dissipation range to other studies from both observations 
and simulations. In the case of the dissipation range, there 
are only really two of these: the study of Al exandrova et al.| 
(2008 ) who used norrnal-mod e (25 Hz) Cluster observations, 
and (Cho & Lazarian| ( 2009 1 who conducted 5 12-' Electron 
MHD (EMHD) simulations. Although EMHD is a high fre- 
quency model for magnetic field dynamics with inertialess 
ions (the opposite of Hall MHD), within the scales that we are 
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analysing here, pj^ '^k< p~^, the induction equation for both 
Hall MHD and EMHD are near-identical for an or der unity 
/3, plasma. Thus the authors ( Cho & Lazarian|200 9') also use 
EMHD to study the dynamics and statistics or magnetic field 
fluctuations in solar w ind turbulence. 

,Alexandrova et al.|(200 8 ) used a continuous wavelet trans- 
form using the Morlet wavelet to calculate both PDFs and flat- 
ness (kurtosis) functions. They too found non-Gaussian PDFs 
within the dissipation range but, in contrast to our monoscal- 
ing results, showed that the dissipation range has a kurtosis 
which increases rapidly with scale r i.e. suggestive of mul- 
tifractal scaling. They also showed a dependence of this flat- 
ness function with /?, . The increase of compressible fluctua- 
tions coincided with a steepening of the spectral slope and the 
increase of intermittency, as reflected in the flatness function. 
This prompted the authors to suggest that, another non-linear 
energy cascade, akin to the inertial range, takes place in the 
dissipation range. However, due to the increased role of com- 
pressible fluctuations at these scales, driven by plasma density 
fluctuations, this is a small-scale compressible cascade. Due 
to the domin ance of the Hall term at these scales, [Alexandrova| 
|et al.| ( [2008| then suggest a simple model based on compress- 
ible Hall MHD, which takes into account density fluctuations 
to describe the deviation of the energy spectrum slope from 
the -7/3 predicted pu rely by the i nduction equation for in- 
compressible EMHD (Biskamp et al. 1996). Without a fur- 
ther study of more data intervals, it is difficult to say why our 
interval shows quite a different scaling behaviour to the in- 
terval investigated by |Al£X£ndrovae^l. (2008). However, 
considering our results and our earlier arguments regarding 
the magnetic compressibility and the role of the Hall term, 
the notion of a compressible cascade is an appealing one. Al- 
though specific results we re not explicitly shown, we should 
also mention the work of Sahr aoui & Gol dstein ("2 010") ) , where 
the monoscaling result presented in this article and in ( |Kiyani| 
|et al. 2009b) was also independently confirmed with Cluster 
magnetic field data. 



The very informative article by Cho & Lazarian ( |2009 1 de- 
scribes the results of their EMHD turbulence simulations, as 
well as ERMHD simulations. The UDWT digital filters used 
in our study are similar (in spiri t) to the multi -point structure 
functions used in the analysis of Cho & Lazarian (2009 ). The 
main feature of these simulations which is relevant to our re- 
sults is that, for the EMHD simulations, the monoscaling be- 
haviour found in our results is also seen in the higher-order 
statistics computed from multi-point structure functions for 
moment orders up to m = 5. Computing higher order moments 
greater than 5 becomes very difficult in non-Gaussian heavy- 



tailed statistic s due to the poor sampling in the tails ( jDudok 
De Wit|2004| ). The scaling exponent H for their <^(m) = Hm 
scaling (as they used C(p)/CO)) is not explicitly stated in 
their paper However, this is easily obtained from the -7/3 
spectral slope that they obtain in their energy spectrum (for 
both EMHD and ERMHD), using the well known relation- 
ship P = -(2H+ 1) ( IKiyani et al.|2009b[ ), where /3 is the spec- 
tral slope. We can then infer that H = 2/3 ~ 0. 66. This is the 
scaUng for rj^ (kj^) variations. The results of |Cho & Lazar- 



|ian| ( |2b09| ) also show moderately non-Gaussian PDFs. The 
authors then compute two-point flatness functions, as well 
as wa velet flatness functions u sing the Morlet wavelet sim- 
ilar to Alexandrova et al. ( 2008) ). Interestingly, for rj^ they 
find a nominal dependence of the flatness with separation 
scale r (equivalent to r) for the two-point flatness function; 
but a larger increase with r for the wavelet flatness function. 



This last part offers slightly conflicting results when consid- 
ering their multipoint structure function monoscaling for mo- 
ments o f order less than 5. Nevertheless, the results of |Cho| 
|& Laz arian ( 2009 ) offer the closest agreement to our scaling 
results presented in this article, and serve as a good compari- 
son to an appropriate model for dissipation range fluctuations 
(EMHD induction equation as opposed to ideal MHD). Sim- 
ilarly to Alexandrova et al. ( 2008) , we could speculate that 
the difference m our exponents, H± = 0.98 and //|| = 0.85 
(and also the respective spectral exponents), and the -7/3 
from EMHD could be due to the differing role of density 
fluctuations. However, in recent gyrokinetic simulations of 



anisotropic plasma turbulence, Howes et al.| ( 2011b) show that 
steeper than -7/3 spectra are only obtained if the full kinetic 
physics is considered in terms of retaining a physical damping 
mechanism e.g. collisionless Landau damping, transit time 
damping etc. This is absent in the EMHD, ERMHD and Hall 
MHD descriptions (see also (Howes 2009)). The gyrokinetic 



descrip tion retains this missing physics and thus [Howes et al. 
( 201 lb) ) claim that their simulat ions succ essfully re cover the 



steeper scaling seen here and in (|Kiyani e t al. 2009b Alexan- 
jdrova et al.|2 009'; Sahraou i et al.|2010F and ( jChen et al.|2010) ) 
for similar plasma parameters. 

3.2.2. Data uncertainties and errors 

There is a small but important data caveat that one should 
mention with regards to these results. This concerns the noise 
floor of the STAFF search-coil magnetometer and the signal- 
to-noise ratio at small scales near the electron gyroradius. Due 
to the dyadic frequency spacing of the UDWT that we use, fre- 
quency resolution is poor. We pay this expense for the faster 
algorithms and the smoother spectra that we gain; the latter 
being preferred for broadband spectra where robust estimates 
of the spectral slope are needed. This gives the impression 
that the signal to noise ratio at the smaller scales is higher than 
it might actually be. In reality there exist some high power 
spikes in the noise-floor signal > 60Hz which could contam- 
ina te the observations at the electron gyro-scale (see figure 1 
in ( Kiyani et al.|2 009b )). These spikes are due to interference 
from electrical signals of the Digital Wave Processor (DWP) 
instrument onboard Cluster (Comilleau-Wehrlin 2012). The 
DWP signals are used to synchronise instrument sampling 
within the experiments of the Cluster Wave Experiment Con- 
sortium (WEC), which STAFF belongs to. These spikes only 
affect the last data point in figures |T| and J2] and thus should 
not change our results and conclusions significantly. Inter- 
vals with higher power such as those in the slow solar wind 
stream would be needed for better estimates at these higher 
frequencies. In addition, the noise-floor that we use is not 
the actual noise-floor but a proxy for it, obtained from a very 
quiet period in the magnetotail lobes where the magnetic field 
power is very low. The actual noise-floor could be smaller A 
discussion on the noise-response of the STAFF search-coil in- 
strument an d its sensitivity in diff erent plasma conditions can 
be found in| Sahraoui et al.| ( |201 1) ). 

We would also like to note that the PDF for the parallel 
fluctuations seen in fig. [4] for the inertial range from ACE, 
show an unusual skew towards positive fluctuations. This is 
an artifact for this particular interval of ACE magnetic field 
data due to large spike-drops (discontinuities) in the magnetic 
field magnitude |B|. As these spikes arise in |B| and not in the 
components, this reflects itself more in the parallel component 
of the fluctuations rather than the transverse components as, in 
many cases, one can use |B | as a proxy for parallel fluctuations 
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e.g. in Alexandrova et al.] ( 2008[ l. This is also reflected in the 
large errors for the parallel scaling exponents in the inertial 
range shown in fig. [3] (a ii.). However, as these spikes are 
primarily in the large rare events in the positive tail of the PDF, 
we believe that the differences in the PDFs of the transverse 
and parallel fluctuations within the inertial range still occur 
and our discussion above remains valid. 

3.2.3. Discussion 

The isotropy that we observe between the power in the dif- 
ferent fluctuation components at kp^ ^ 1 is by no means a 
result which would apply to all plasma environments. In- 
deed, there is some strong observational evidence that the 
ion/proton plasma /? ,■ also has a strong rol e to play here. 
This was shown by [Hamilton et al.' ("2008^, who repeated 
the inertial range analysis of |Smith et al. (2006b ) within the 
dissipation range and showed that the variance anisotropy 
followed the empirical power-law relationship 6B^/SbJ^ = 
PSD^/PSD\\ ^ /3p"^^ in open field line intervals of the so- 
lar wind from the ACE spacecraft. However, these results 
were limited to 3Hz cadence and thus could not sa mple as 
much of t he dissipation range as in this article. Also Harml-| 
|ton et al.] pOOS ) used a global field (mean field over the en- 
tire mterval) rather than the local scale-dependent field used 
in this article. This same dependence of the magnetic com- 
pressibility on /3, was also indica ted in higher cadence (25 
Hz) Cluster data in the work by Alexandrova et al!] ( 2008 1, 
again suggesting that higher /3, implies higher levels of Cy 
(these authors used (5|B| as a proxy for 5B\\). The effect of 
varying /3, is not mentioned explicitly in our arguments us- 
ing the Hall term explained above, as we have neglected den- 
sity fluctuations to simplify the discussion. However, purely 
from the induction equation point of view, /3, is taken into ac- 
count in dissipation range scales if we retain the effects aris- 
ing from density fluctuations at scales kp, > 1 as in the for- 
mulation ofElecto^ (Schekochi- 
|hin et al.||2009) |Cho & Lazanan||200 9l). In ERMHD the ex- 
pUcit /3, dependence is shown, assuming pressure balance, 
in the pre-factor of the Vj^ x (B • VBj^) term in the induc- 
tion equation for B\\. Thus, this addition of density fluctu- 
ations can generalise our arguments above which were es- 
sentially based on the Hall term alone. At higher or lower 
than unity /3,, this will break the isotropy of the fluctuation 
components seen in this article at kpe ^ 1, as an additional 
anisotropic source of magnetic compressibili ty will be intro- 
duced via the density fluctuations (see also (Malaspina et aLJ 
|2010| l for high cadence measurements of density fluctuations 
in the solar wind between pj^ < k < p~^). The explicit /?,■ de- 
pendance for the rising compressibility is also seen in the fluid 
version of the kinetic Alfven wave solution of compressible 
Hall MHD where the compressible component to the linean 
solution goes as (5B|| ^ A:A,(/3, + 1 + sm^ 9^8)/ sinOtB (Sahraoui 
|2003|, where A, is the ion inertial length, and O/tB is the an- 
gle between the wavevector k and the magnetic field vector 
B. This solution assumes that sind/tB 7^0; thus it is not valid 
for purely parallel propagating Alfven waves. However, the 
rising compressibility is clearly seen here as A: — > A,-. 

Note that in most of the previous arguments we have been 
slightly cavalier in our discussion of field components. Our 
introduction of the local background magnetic field should 
not be looked upon lightly. It is clear that in cases where 
6B <C Bo, where Bo is the global background magnetic field, 
both global and local field descriptions should have negligi- 
ble differences. However, as we observed in the section de- 



scribing higher-order statistics, it is not obvious that this is 
the case here, as the distribution of SB is highly non-Gaussian 
with very heavy tails. In our explanation of the rising com- 
pressibility in sub-Larmor scales using the Hall-term and also 
the lower levels of the compressibility in the inertial range us- 
ing the convective term in the induction equation, we did not 
mention the spatial (and possibly temporal) gradients which 
arise due to the now locally varying background magnetic 
field. In a global DC background field, or very slowly evolv- 
ing field, these gradients are null or negligible. If there is a 
large spectral separation between the scales of the background 
magnetic field and the fluctuations that we are projecting on 
it, then this conclusion could still be quite valid. However, as 
our UDWT method assures that no such large gap exists, the 
statistics we calculate will include the effects of such gradi- 
ents. Further investigation and mathematical grounding of a 
local-field is of crucial importance, as the notion of using a 
local versus a global background magnetic field is intimately 
tied with the issue of whether we have a strong background 
guide fi eld or a weaker more stoc hastic one (see the discus- 
sion by |Schekochihin et al. ( |2009 1 and references therein for 
further details and exploration of this topic ; and also the re- 
cent discussion by Matthaeus et al.| ( [20T2 |) on the stochastic 
nature of the local background field and its relationship with 
global background field statistics). This also brings into ques- 
tion the use of linear wave modes, such as the kinetic Alfven 
wave, which are normally described as small amplitude per- 
turbations on a static or slowly evolving equilibrium field. Al- 
though finite amplitude propagating nonlinear Alfven waves 
are an exact (Elsasser) solution to the ideal MHD equations; 
and similar arguments in the c ontext of the Electron Red uced 
MHD equations are made by Schekochihin et aL] ( 2009) for 
finite amplitude kinetic Alfven waves; this topic also requires 
further exploration (see also the excellent discussion of this 
by Dmitruk & Matthaeus (2009), who explore the importance 
of the strength otHie mean background field using direct nu- 
merical simulations of the incompressible MHD equations). 

One might question a conclusion of a single physical pro- 
cess or type of mediating dynamics (in the dissipation range) 
by the observation that both transverse and parallel compo- 
nents have slightly different scaling exponents and spectral 
indices. This can be reconciled by realising that there is no 
contradiction in the fact that the two components can have 
slightly different scaling exponents but have identical PDF 
functional forms. Take for example the case of fractional 
Brownian motion (fBm) which is a non-Markovian general- 
isation of Brownian Motion. Two fBm stochastic processes 
with different scaling exponents will still have the same Gaus- 
sian functional form for their fluctuations. In fBm the scal- 
ing exponents are simply an additional parameter representing 
the degree of statistical depen dence betw een the fluctuations 
(jMandelbrot 1983; Samorodnitsky^ Taqqu 1994). 

If we look from the point of view of linear wave mode char- 
acteristics, the distinctness between inertial and dissipation 
range scalings could also be seen to reflect the decoupled and 
coupled nature of the wave modes in the inertial r ange and dis- 
sipation ranges respectively ( [Leamon et al. 1999| ). The decou- 
pling of these modes in the inertial range also results in the de- 
coupling bet ween transverse and compr essible magnetic field 
fluctuations ( [Schekochihin et al.||2009] ). In this picture the 
scaling in the dissipation range could then be described by 
non-linearly interacting kinetic Alfven waves. In all these 
cases, as they involve linear wave modes, the non-linear cas- 
cade would then take the form of 'critical-balance' type cas- 
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cades as in the work s of |Schekochihin et al.| ( |2009l l and |Howes| 
et al. ( 2008 2011b| l; or if one can show that the non-hnearity 



IS sm all, a weak-turbulence' type description (Galtier et al. 
|2002| l. Although there is no justification for the latter in the 
inertial range, there might be a case for it in the dissipation 
range if the coupling between modes becom es weak and re - 
sults in small amplitude fluctuations (Ruda kov et al.|[T011| l. 
In hydrodynamics described by the Navier-Stokes equations, 
this latter assu mption can be va lid if one uses a local Reynolds 
number (see ( |Batchelor||1953[ )), which is small in the high-A: 
regime (when the spectral amplitudes become small) well into 
the hydrodynamic dissipation range. This can then provide a 
suitably small parameter to base a convergent (or asymptotic) 
perturbation expansion. Of course, without concrete mea- 
surements of non-linearity, this is all speculation. However, 
it could possibly be an avenue for a future (maybe fruitful) 
investigation. 

An open question is why such a transition from multifractal 
scaling to self-similar mon scaling occurs between the iner- 
tial and dissipation ranges. Kiyani et al. ( 2010 1 briefly dis- 
cuss this using hydrodynamic analogies (Frisch & Vergas- 



|sola|199l] Chevillard et al. 1200 5) of a near-dissipation rangi 
where the increased effect of viscosity successively switches 
off the available scaling exponents of the multifractal field - 
equivalent to narrowing the '/(a)' spectrum in multifractal 
models . However, no such apparent analogies with a vis- 
cous 'fractal-dampening' can happen in a near-collisionless 
solar wind. Also, unlike the hydrodynamic analogy this tran- 
sition happens very fast as we cross kpt ^ 1 ; although a scal- 
ing analysis within the 'transition range' might yet reveal an 
intermediate scaling behaviour which shows this successive 
'switching-off of scaling exponents. Indeed, the steepening 
of the spectral slope in the transition region is strongly sug- 
gestive that such an intermediate scaling range could possibly 
exist. This is an exciting prospect and will be the subject of 
a future investigation. This latter investigation is not feasible 
with our UDWT technique as due to the dyadic spacing, it 
has poor frequency resolution - a more detailed 'fine-toothed 
comb' technique such as the use of wavelet packets (WaldenI 
& Crista n 1998) or the continuous wavelet transform (Hor-| 



bury et al.|2008t|Podesta.j2009| would be needed here. 



4. CONCLUSIONS 

In this article we have performed a systematic scale-by- 
scale anisotropic decomposition of magnetic field fluctua- 
tions in the solar wind across both the inertial and dissipation 
ranges, from MHD to electron scales. We use a novel scheme 
based on discrete wavelet filters which self-consistently de- 
composes the magnetic fields into scale-dependent back- 
ground and fluctuating fields. From the background fields 
we then construct a local orthonormal scale-depedent field- 
velocity coordinate system, on which we then project the fluc- 
tuations and calculate our statistical quantities. In summary, 
our main findings are:- 

1. There exists a successive scale-invariant reduction in 
the power ratio between parallel and transverse com- 
ponents as we move to smaller scales below the ion- 
Larmor radius p,-. Importantly, at kp,, ~ 1 with k 
components highly oblique/perpendicular to the back- 
ground magnetic field, the fluctuations become isotrop- 
ically distributed between the power in all three compo- 
nents (two transverse and one parallel) i.e. equipartition 
of magnetic energy between all components. 



2. By comparing with linear solutions of the Vlasov equa- 
tion and computing the magnetic compressibility we 
show that the scale-invariant reduction in the power ra- 
tio is also qualitatively consistent with the transition to 
kinetic Alfven wave solutions which rapidly develop 
magnetic compressible components at scales around p, . 

3. Using a Hall MHD model we provide a simple (non- 
linear) description of this new result in terms of how 
the Hall term is responsible for the rise of parallel 
(compressible) magnetic fluctuations and the resultant 
isotropy at scales around p^,. 

4. Higher order statistics and PDFs for parallel/transverse 
fluctuations reveal that the dissipation range is char- 
acterised by a single monoscaling statistical signature 
with an isotropic distribution of fluctuations; thus sup- 
porting the above observations of isotropy. This is in 
contrast to the inertial range where the corresponding 
signature is anisotropic and multiscaling (multifractal). 

In the compressible Hall-MHD framework, at such scales 
the velocity field is expected to decouple from the evolution of 
the magnetic field and some researchers ( Servidio et al. 2007J 
Alexandrova et al. 2008j ) have suggested that a new non-linear 
compressible cascade is responsible for the new power-law 
regime in the PSD. In this context, we should also mention 
that all of the arguments presented in this paper do not in- 
clude any discussion of the possibility of phe nomena arising 
from the effects of (turbulent) reconnection ( |Lazarian et al | 
2012| l or of the possible role of a turbulent or renormalized 
magnetic diffusivity - of which either, or both, are purported 
to be important at such scales. Neither have we discussed the 
role of density fluctuations, which have been neglected in our 
arguments, but are shown to play a p ossibly significant role 
through the electron compressibility ( [Gary & Smith 2009). 
Future work would involve looking at the possible importance 
of such fluctuations as well as investigating the energy trans- 
fer between modes in Hall-MHD/Electron MHD theories. 

In the context of theories of turbulence which invoke 
plasma wave modes as the mediators of the energy cascade, 
two candidate wave modes are classically suggested in the 
dissipation range: quasi-perpendicular kinetic Alfven waves 



and quasi-parallel Whistler waves. As first suggested by Gary 



& Smith (2009) and shown in I Salem et al. (2012), the mag- 
netic compressibility, C||, might be a more robust measure 
ment to distinguish between the different wave modes possi 
bly mediating the energy c ascade, than the phase speed plots 



calculated hy E/B ratios (Bale et al. 2005 [ Sahraoui et al. 
2009 1. In hot plasmas (/3, ^ 2) and at highly oblique angles 



to the background magnetic field, classical Whistlers arising 
as an extension of the fast magnetosonic modes at high fre- 
quencies > LJci (where ajd is the ion cyclotron frequency), 
are strongly damped and s plit up into different io n Bernstein 
modes. Recently, however, Sahraoui et al. ( 2012| l have shown 
that, under the same condition, a Whistler-like branch hav- 
ing similar properties (e.g. dispersion and polarisation) as the 
classical Whistler mode, rises from the kinetic Alfven wave 
mode and extends it to very high frequencies cj ^ w„. Also, if 
there is a significant amount of power in shallow angles, e.g. 
30°, between k- vectors and the background magnetic field, 
then the classical quasi-parallel Whistlers could also show the 
same amount of magnetic compressibilit y that we observe in 
this paper ( .Gary & Smith 2009, Salem et al.|2012| l. 
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Our study of higher order statistics show a new and sig- 
nificant result of the change of anisotropy at p,. As in the 
behaviour of the inertial range, these scalings and statistics 
are a distinctive and important characteristic of the turbulence 
at sub-ion Larmor scales and any cascades which take place 
there. Therefore it is essential that any comprehensive statis- 
tical theory of plasma turbulence for such scales, be able to 
describe such statistics. 

An open question is whether these results depend upon bulk 
plasma parameters such as /3, and, in particular on the angle 
of k to the background magnetic field. A limitation of the cur- 
rent work is that it is based on a handful of data sets with very 
similar plasma parameters. A significant and maybe conclu- 
sive advance would be made if a more comprehensive survey 
was conducted. This would undoubtedly be a fruitful avenue 
for future investigations. 
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APPENDIX 

A. ORTHONORMAL SCALE-DEPENDENT BASIS 
AND THE UNDECIMATED DISCRETE 
WAVELET TRANSFORM 

We use the Undecimated Discrete Wavelet Transform 
(UDWT) method (also known as the stationary, tr anslation- 
invariant, r edundant or a 'trous wavelet transform ( |Nason &| 
|Silverman| pr995; Ogden 1997, Mallat 2009)) to decompose 
the vector magnetic field observations into scale dependent 
background and fluctuating components. Unlike the stan- 
dard discrete wavelet transform (DWT), the UDWT does not 
downsample the data at each stage of the transform. Thus, the 
UDWT provides information about the signal at e ach observa- 
tion time a nd so reta ins event information ( W alden & Cristan| 
[T998, Percival & Walden 2000). This timing information is 
necessary to project the fluctuations onto a mean background 
field at each observation time. However, the extra redundancy 
due to no downsampling means that the orthogonality proper- 
ties of the DWT are no longer present in the UDWT. 

In this article we use the Coiflet 2 wavelet (coifZ), which 
provides a 12-tap high/low pass filter pair, chosen due to the 
compromise between time-frequency compactness, smoother 
spectral index estimate, and as its first fou r moments are zero 
it captures spectral indices as steep as -9 (|Farge & Schneider 



[2006 ). This latter property is very important as the standard 
method of using increments (equivalent to using a Daubechies 
order 1 (dbl), or Harr, wavelet with only one zero moment) 
to calculate two-point structure functions in turbulence is lim- 



ing the mertial range o: 



ited to spectral slopes shallow er than -3 ( [Farge & Schneider 
2006] |Cho & LazaiTan 
the mei 



2009 1. This is adequate for study- 



> k approximations 
(large scales) 



a1 



details 
(small scales) 

d1 

¥- 




frequency 



Figure Al. (Upper panel) A schematic, in terms of spectral frames, showing 
the first stage of a UDWT decomposition in frequency space: the raw signal 
is decomposed into a low-pass (approximation), 'al', and a high-pass (detail) 
component, 'dl', using the coifi wavelet filters (low-high pass quadrature 
mirror filter pair). Each of these components has an equal (half) share of the 
frequency space (Lower panel) At the next stage, the 'al' low-pass filtered 
signal is then further decomposed, into a low-pass signal, 'a2', and a high 
pass-signal, 'd2' (again equal share of frequency space). In this way the sig- 
nal is decomposed into dyadic frequency bands - equivalent to the application 
of a dyadic filter bank of band-pass filters to the signal. The resultant set of 
'details' {dl,d2,d3 ...) provide the fluctuations, and the set of 'approxima- 
tions' {al,a2,a3 ...) (once inverse transformed) provide the respective back- 
ground fields to be projected upon. The stars represent the central frequency 
of the resultant band-pass filtered signals used for the fluctuations. 



is ubiquitously less than -2 ( [Salem e t al."2009), but is inade- 
quate in the dissipation range where spectral slopes have been 
observed (Leamon et al. 1998ai Smith et al. 2006a) and the- 
oretically predicted ( Schekochihin et al.^2009) to have much 
steeper slopes between -3 and -5. Crucially, for data analy- 
sis purposes, the filters corresponding to the coif2 wavelet are 
virtually symmetric and thus possess a virtually linear phase 
response (Daubechies 1992). The coif2 wavelet (and corre- 
sponding scaling fun ction) filters used here are identical to 
the ones published inlDaubechies ( 1992). 



The UDWT scheme presentedoel ow is sim ilar to the 
wavelet technique pioneered by Horbur y~et al.| (f2008 ) who 
used a continuous wavelet transform (CWT), where the back- 
ground magnetic field was defined as convolutions with the 
scale-dependent envelope of the Morlet wavelet used. How- 
ever, the method described here differs from the CWT scheme 



turbulence where the spectral slope 



used by [Horbury etaL| ( |2()08| ), as the DWT and UDWT filters 
used are explicitly designed self-consistently to provide both 
fluctuations and background fields such that no information is 
lost and the signal can b e reconstructed exa ctly - a so-called 
quadrature mirror filter ( [Daubechies || 1 992|l . Thus, these fil- 
ters also ensure that there is no spectral gap between the back- 
ground field and fluctuations. Details of the actual fast pyra- 
mi dal algorithms used to im plement the UDW T can be found 
in ( [Percival & Walden|200 0) and (Manat 2009). For this arti- 
cle we simply attempt to summarise what are the salient out- 
comes of the technique. Thus, the description will be neces- 
sarily brief with limited details. A detailed description of the 
technique as applied to plasma turbulence anisotropy, as well 
as an extension of this to non-dyadic wavelet packets, will be 
described in a longer manuscript. 
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The UDWT is equivalent to the following operation: at each 
temporal scale r (corresponding to a frequency) we succes- 
sively apply a pair of high and low -pass wavelet filters on the 
data: these provide the 'details' (fluctuations) and the 'ap- 
proximations' (precursor to background field). Since these 
are dyadic filters, they divide the available frequency space 
in half at each stage. These filters are then upscaled by two 
and the process is then repeated on the previous approxima- 
tion signal. The results of the first two stages of this process 
are shown in fig |Al| Importantly, at each stage of the UDWT, 
it is the approximations which are being decomposed into a 
further set of approximations and details. The effective end 
result of this successive filtering is a series of dyadic band- 
pass filtered signals of the magnetic field R(tj) (sampled at 
discrete times tj) to give a fluctuation in terms of the wavelet 
coefficient 6B(tj,T) (at time tj based on a temporal scale r): 



6B(tj,T) = ^B(r,) [V^^(Th-tj)] 



(Al) 



k=l 



where r = 2'A : / = {0, 1,2,3, . . .} is the dyadic scale param- 
eter, A is the sampling period of the observations, and V' is 
the effective band-pass filter comprised of a succession of the 
low-pass and high-pass filters. Note that without the factor 
^/T and choosing th e db l (Harr) wavelet instead of the coif2, 
the definition in eq. ( |A1| | is identical to calculating increments 
which are used in structure function analysis of intermittency 

^1 



studies in turbulence (Katul et al. 2001 
[20061 ISalem et al.|2009| l. Note also that 



Farge & Schneider 
e resultant wavelet 



coefficients dR(tj,T) have to be phase-corr ected due to the lin- 
ear ph ase shift introduced by the filtering ( [Walden & Cristan] 
|1998| l. Also, edge artifacts due to incomplete convolutions of 
the signal with the filters, are discarded before the construc- 
tion of the 'approximation' and 'detail' signals. 

The scale r can be related to a central frequency / in Hz 
( Abry|1997| of a dyadic frequency band and so we can write 
oB{tj, f). The result of the low-pass filter at each stage can 
similarly be written as B(f,,/) using the wavelet conjugate 
scaling function (Ogden,^1997). It is important to note that 
both SB(tj,f) and B(fy , /) Uve in a function space spanned by 
the bases constructed from shifts and dilations of the coif2 
wav elet a nd scaling functions (i.e. shift tj and dilation r as in 
eq. (Al i). We use the inverse UDWT operation on B(fy,/) 
in order to obtain the mean background field B{tj, f), in the 
direction of which we then project the fluctuations SB(tj,f) . 

The mean field direction unit vector is then defined as 
e||(/y,/) = B(tj,f)/ |B(fy,/)|. This is distinct from the mean 
field coordinate pre sented m (Belcher & Davis 1971 Lea- 
|mon et "aL] |1998b| l as the mean field here is a locally 
varying scale-dependent field consistent with the scale de- 
pendent fluctuations. The other two perpendicular axes 
are ti^yitjj) = (e||(f;,/) x (Vsw)) / |e||(f/ J) x (Vsw)| and 
e_L2(f;)/) = e||(0,/) X e_Li(fj, A where (Vsw) is the solar 
wind velocity direction unit vector time-averaged over the en- 
tire interval. Since the solar wind velocity is in a very fast 
and steady stream, and is within ^ 3° of the GSE x direc- 
tion in both the ACE and Cluster intervals, it is reasonable 
to take a time averaged global, rather than a local, velocity 
field. Together [e±i(tj,f),e±2itj,f),e\^(tj,f)] form a time and 
scale dependent orthonormal basis. The fluctuations paral- 
lel 5B]^\{tj,f) and transverse [i5B_Li(/j,/) , (5B_L2(fjj/)] to the 



mean field are then given by projections onto this new basis. 
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